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ABSTRACT 

We present the complete galaxy cluster catalog from the Northern Sky Optical Cluster Survey, a 
new, objectively defined catalog of candidate galaxy clusters at z < 0.25 drawn from the Digitized Sec- 
ond Palomar Observatory Sky Survey (DPOSS). The data presented here cover the Southern Galactic 
Cap, as well as the less-well calibrated regions of the Northern Galactic Cap. In addition, due to 
improvements in our cluster finder and measurement methods, we provide an updated catalog for the 
well-calibrated Northern Galactic Cap region previously published in Paper II. The complete survey 
covers 11,411 square degrees, with over 15,000 candidate clusters. We discuss improved photometric 
rcdshifts, richnesses and optical luminosities which are provided for each cluster. A variety of sub- 
structure measures are computed for a subset of over 11,000 clusters. We also discuss the derivation 
of dynamical radii r2oo and its relation to cluster richness. A number of consistency checks between 
the three areas of the survey are also presented, demonstrating the homogeneity of the catalog over 
disjoint sky areas. We perform extensive comparisons to existing optically and X-ray selected cluster 
catalogs, and derive new X-ray luminosities and temperatures for a subset of our clusters. We find that 
the optical and X-ray luminosities are well correlated, even using relatively shallow ROSAT All Sky 
Survey and DPOSS data. This survey provides a good comparison sample to the MaxBCG catalog 
based on Sloan Digital Sky Survey Data, and complements that survey at low redshifts 0.07 < z < 0.1. 
Subject headings: catalogues - surveys - galaxies: clusters: general - large-scale structure of the 
Universe 



1. INTRODUCTION 

The construction of large catalogs of galaxy clusters 
for use in studies of cosmology, large scale structure, and 
gala xy evoluti on has often proven to be a difficult task 
(see lGa]ll2(M for a review). Indeed, the last such cat- 
alog generated using optical data, and covering t he en- 
tire hi gh-galactic-latitud e Northern sky was that oflAbel 
([l958h . updated in 1989 (lAbell. Corwin fc Qlowini ri989L 
More recently, the Nort hern ROSAT All-Sky (NORAS) 
Galaxy Cluster Survey l|B5hringer et al.1 [20001 ) has pro- 
vided an X-ray selected catalog covering a similar region, 
but with many fewer clusters, wh ile the largest catalog 
using purely digital observations (jKoester et al.l l2007ri ) 
covers a smaller area, albeit with much better photome- 
try. 

Because much improved optical data has become avail- 
able, with automated techniques to generate objective, 
well-characterized cluster samples, we undertook the gen- 
eration of a modern, optically selected cluster cata- 
log, the Northe r n Sky Optical Cluster Survey (NoSOCS 
iGal et all l2000l I2003L hereafter Paper s I and II), and 
its deeper extension ( Lopes et al.ll2004l hereafter Paper 
IV). The need for a modern cluster survey covering a 
significant portion of the sky is striking. The catalog 
from Paper II has already been used to suggest a connec- 
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tion between short-durat i on gamma-ray burs ts and clus- 
ters dGehrels et all 120051: iBloom et al.l [2006h . to search 
for giant arcs (He nnawi et a l. 2008), to associate com- 
pact groups and large-scale structure (|de Carvalho et all 
I2005t lAndernach fc Coziolll2006ft. and to ex amine X-ray 
and optical cluster properties (jLopes et al.ll2006| ). Here 
we present the second and final installment of this cat- 
alog, including photometric redshifts, richnesses, optical 
luminosities and substructure measures. Although supe- 
rior imaging data is now available from the Sloan Digital 
Sky Survey (SPSS. lYoTk et al.ll2000T) . the imaged area 
covers less than half of the Northern sky, and is not 
expected to ever reach the area coverage of our cata- 
log. The currently pu blished cluster catalog from SDSS 
(|Koester et al.l l2007bf ) covers ~ 7500 deg 2 , containing 
nearly 14,000 clusters. That survey, using the MaxBCG 
technique, has a lower redshift cutoff of z = 0.1; our sur- 
vey extends down to z = 0.07. This provides a sample of 
more l ocal clusters whose properties can be examined in 
detail (jLopes e t al. 2008), especially using the extensive 
SDSS spectroscopic database. 

The regions covered in this paper are shown in Fig- 
ure [1] Three separate areas are covered: (1) The well- 
calibrated North Galactic Pole (NGP) region already de- 
scribed in Paper II (dotted lines); (2), the portion of 
the NGP with less well calibrated plates, not covered 
in Paper II (solid lines); and (3), the Southern Galactic 
Pole region (SGP, dashed lines). The area covered by 
the NGP-poor region is 2813D°, while the SGP covers 
2917D . Together with the 5681D° surveyed in Paper II, 
the final NoSOCS catalog covers 11411 square degrees. 
The distribution of clusters in the survey is shown in 
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Figures [2] and [3l in equatorial and galactic coordinates, 
respectively. Only clusters with richness N ga i s > 20 arc 
shown for clarity. 

The survey methodology is described in §2. Although 
the overall detection technique is very similar to that dis- 
cussed in Paper II, we have modified our definition of bad 
areas due to very bright objects, reducing contamination 
by spurious detections. Changes to the photometric red- 
shift estimation yield more robust (and realistic) errors. 
The richness estimator has also been improved, provid- 
ing robust error estimates, as well as total r-band optical 
luminosities. Then, in §3, we discuss two new sets of pa- 
rameters computed for our cluster sample: estimates of 
cluster substructure and the dynamical radius r2oo- 

In §4 we describe the general characteristics of our clus- 
ter sample, and present consistency tests for the three 
sky regions utilized. The selection functions describing 
the completeness as a function of richness and redshift 
are presented, as is an estimate of the contamination by 
projection effects. Our complete and final cluster cata- 
log, including an updated version covering the area from 
Paper II, is presented. This cluster sample is compared 
to other optically-selected catalogs in §5. In particular, 
w e compare our catalog to the SDSS MaxBCG catalog 
of iKoester et al.l (|2007bf ). the only modern optical clus- 
ter catalog covering a similar area and redshift range. 
We then examine the correlation between our redshifts 
and richnesses and the X-ray measurements from the 
North ern ROSAT All-Sky G alaxy Cluster Survey (NO- 
RAS, [BShringereEj^EOOO). Because NORAS consists 
of many fewer clusters than our catalog, we use our op- 
tical positions to measure X-ray fluxes and upper limits 
from RASS for a significant subsample of NoSOCS. 

2. SURVEY METHODOLOGY 

The detection of galaxy clusters in modern optical 
imaging surveys typically utilizes the existence of a 
tight color-magnitude relatio n for cluster g a laxies, noted 
nearl y half a century ago (|Bauml Il959t iBower et al.1 
1992). Surveys at low and high redshift have made 
use of this observation to efficiently detect clusters, 
while achieving low contamination (false p ositive) rates 
dGladders k Yedl2000HHansen et al.ll2005h . We use the 
galaxy catalogs from the Digitiz ed Second Palomar Ob- 
servatory Sky Survey fDPOSS. [Diorgovski et"aTlll999h 
as the basis for our survey. Unfortunately, the lim- 
ited photometr ic accuracy (a mag ~ 0.25 m at r = 19.5, 
iGal et al.ll200l of DPOSS forces us to rely solely on the 
two-dimensional projected galaxy distribution for clus- 
ter detection. Details of the photometr ic calibration and 
star /galaxy separation are discussed in lGal et al.l (|2004l ) 
and lOdewahn et all (|2004h respectively. In brief, the po- 
sitions of the galaxies ar e used to genera te adaptive ker- 
nel (AK) density maps (|Silvermanlll986l ) which outputs 
images in u nits of projected galaxy d ensity. We then run 
SExtractor (|Bertin fc Arnoutall996[ ) on these images, de- 
tecting peaks which are identified as potential galaxy 
clusters. We refer the reader to Papers I and II for more 
comprehensive descriptions of the cluster detection. Pho- 
tometric redshifts are estimated using the background- 
corrected mean r magnitude and median g — r color of 
the galaxies within a 0.5ft.~ 1 Mpc radius of the cluster 
center; details on improvements to this estimator from 
our previous work are described below. The photometric 



redshifts are used to recenter the clusters, also discussed 
later. Richnesses are computed by counting galaxies with 
M* - 1 < M* < M* + 2 within the same radius, with 
corrections applied for higher (and now lower) redshift 
clusters where the faint (or bright) end of this magni- 
tude range is beyond our catalog limits. We note that 
galaxy colors are used only in the post-detection steps to 
estimate photometric redshifts for the clusters. 

2.1. Enhanced Removal of Bad Areas 

In the process of generating the catalog presented in 
Paper II, it became obvious that the DPOSS catalogs 
were insufficient for finding very bright stars and galax- 
ies, which are typically deblended into numerous fainter 
components. In later work, Paper IV utilized the Tycho- 
2 catalog to excise candidate higher-redshift clusters in 
the area of bright stars after detection. To avoid spu- 
rious cluster detections due to these artifacts, we now 
rely on the Tycho-2 and RC3 catalogs to exclude regions 
in the vicinity of bright objects before performing clus- 
ter detection. Specifically, around bright stars we ex- 
clude circular regions whose area depends on the star's 
Tycho magnitude; 2' radius for mxycho < 7.0, 1.5' for 
7.0 < mTycho < 8.0, and 1.0' for 8.0 < mT yc ho < 9.5. 
These radii were chosen by visually inspecting plate im- 
ages and the resulting galaxy catalogs to determine the 
sizes of regions contaminated by the bright stars. In ad- 
dition, l arger regions around bright galaxies in the RC3 
catalog (jde Vaucouleurs et al.lll99lT ) were excised, corre- 
sponding to 5 x r RC 3 for r RC 3 < 25" and 8 x r RC3 for 
mca > 25". 

The removal of objects in these contaminated re- 
gions results in empty holes in the galaxy catalogs. 
An undesired consequence is that in these areas the 
adaptive kernel artificially increases the smoothing ra- 
dius (which is inversely proportional to the local den- 
sity). To remedy this problem, we generate simu- 
lated galaxy catalogs using the Raleigh-Levi distribution 
(|Postman. Lauer. Oegerle. Sz Donahue! l2002t Paper II), 
and use these to fill the excised areas as well as the den- 
sitometry spot region on each plate. With this technique, 
the regions which would otherwise be empty instead con- 
tain galaxies with the average projected density for each 
plate. 

After cluster detection is completed, the final catalog 
is again checked against the list of bad areas. For some 
very extended, nearby galaxies, or very bright stars, we 
found through visual inspection that the exclusion area 
described above was not always sufficient. We therefore 
increased the radius of these areas, by a factor of 1.5 for 
the Tycho-2 stars, and to 10 x r R c3 for the RC3 galax- 
ies. In addition, we created additional bad areas from 
catalogs of Galactic globular and open clusters, as these 
objects were found to be contaminants in our previous 
catalogs. Any candidates found within an exclusion ra- 
dius defined by the sizes of these Galactic clusters was 
also eliminated. Visual inspection of candidates flagged 
in this last step shows that the vast majority were indeed 
bad. A total of 404 candidates, 2.5% of the sample, are 
removed from the catalog in this step. 

2.2. Photometric Redshift Improvements 

A number of small but significant changes were made 
to the algorithm presented in Paper II. These modifi- 
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FlG. 1. — The distribution of DPOSS plates used in NoSOCS, showing the NGP well-calibrated (solid lines), poorly calibrated (dotted 
lines) and SGP (dashed lines) regions. 




cations enhance the photometric redshift measurements 
based on DPOSS photometry and provide more robust 
error estimates. First, we re-derived the empirical rela- 
tion between the median g — r color, mean r magnitude of 
the cluster galaxies, and the spectroscopic redshift (equa- 
tion 1 in Paper II). This was prompted in part by modi- 
fications to the background estimator (described below) , 
which results in changes to the global properties of the 
cluster galaxy populations. Our larger sky coverage also 



allowed us to restrict the spectroscopic cluster sample to 
those clusters with more than three concord ant redshifts 
in the compilation of lStruble &: Roodl (|1991h . resulting in 
a training sample of 254 clusters. We also found that the 
redshift estimator was more reliable when restricted to 
galaxies with m r < 19.5, as opposed to m r < 20.0 used 
before. The recalibrated photometric redshift relation 
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used here is 

z phot = 0.3273 X(g- r) med ~ 0.0702 x r mean - 1.2685 

(1) 

with a Qa(z sp ect - Zphot) I '(1 + z S pec) = 0.023, an improve- 
ment of ~ 30% over the results in Paper II. 

As noted above, one major modification to our tech- 
nique involves measurement of the fore- and background 
galaxy contamination in the cluster area. In Papers I 
and II, we used color and magnitude distributions from 
each DPOSS plate (~ 30 deg 2 ), scaled to the area of 
each cluster on that plate, as the background correction. 
This ignores the contribution of local large scale struc- 
ture to the background of each cluster, which can intro- 
duce systematic errors since galaxy colors and l uminosi- 
ties are strongly correlate d with local density (|Dresslerl 
Il980t iBlanton et al.ll2005f ). We therefore implemented a 
local background estimator, as follows: 

1. A random position is chosen within a background 
annulus of width 1.3° starting Sh^ 1 Mpc from the 
cluster center. 

2. A box of size 20' x 20' is placed at the random 
location 

3. A check is performed to see if the box intersects a 
bad area (hole) in the survey; if so, we return to 
step 1. 

4. The distribution of colors and magnitudes is gen- 
erated for the galaxies in this box. 

5. The procedure is iterated until ten background re- 
gions are successfully measured. 

6. The 3cr-clipped medians of the distributions from 
the ten background regions are used as the back- 
ground correction for that cluster. 

The redshift estimator is run ten times for each cluster 
candidate. Changes in the placement of the randomly 
located background measurement regions result in vari- 
ations of the background galaxy color and magnitude 
distributions, which effect the final photometric redshift. 



By repeating the measurement, we derive an estimate of 
the random error in z p ) lot due to the background cor- 
rection, which we then add in quadrature to the scat- 
ter from the redshift - photometric properties relation. 
Although some of the latter is likely due to difficulties 
with properly estimating the background contribution, 
we prefer to estimate the redshift errors conservatively, 
adding the errors as if they were independent. We note 
that such problems are likely to be significantly reduced 
in modern digital imaging surveys, where the photometric 
errors are an order of magnitude smaller than for D P O S S 
(|Gal et all 120031 lAdelman-McCarthv et al.ll2007t ). Asa 
final change to the redshift estimator, we allow more it- 
erations for convergence (15 instead of 10). This was 
found to increase the number of successful redshift esti- 
mates, while allowing further iteration simply grew the 
computational requirements with little improvement. 

For some clusters, the photometric redshift estimator 
does not converge in all ten of the runs. The final pho- 
tometric redshift is taken as the mean z p hot from the 
n S uccessfui runs. The rcdshifts and their associated errors 
are provided in Columns 5 and 6 of Table [3] For those 
clusters where the z p hot estimate always failed, these 
three columns are all set to 0. Such clusters are likely to 
be spurious detections or have contaminated photometry 
from bright stars, telescope reflections or other artifacts. 

2.3. Cluster Centroids 

During the photometric redshift measurement process, 
the cluster positions are also recomputed. At each iter- 
ation of the Zphot computation, we calculate the median 
position of the galaxies within a lft. -1 Mpc radius of the 
previously determined center, and this is taken as the 
new cluster centroid for the next iteration of the pho- 
tometric redshift estimation. To avoid large offsets, the 
maximum change in position is limited to 2 arcminutes; 
if the recomputed center is further from the previous lo- 
cation, then the center is not moved. The above steps 
are repeated for each of the n success f u i photometric red- 
shifts, and the final cluster position is recorded as the 
mean of the corresponding n success f u i positions. In those 
cases where the photometric redshift estimator does not 
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converge, we retain the original position (from running 
SExtractor on the density map). 

2.4. Richness and Luminosity Measures 

Richnesses and luminositie s are compu t ed us ing the 
basic methods described in lLopes et al.l ((2006), with 
some further refinements, and adjusted for the lower red- 
shift range probed here. The procedure consists of five 
steps for each cluster, described below. 

1. We use Zphot to determine the apparent magnitude 
m*, the aperture corresponding to 0.50 hr 1 Mpc, 
and the /c-corrections k e and k s for elliptical and 
late- type galaxies (Sbc) at the cluster redshift. We 
select all galaxies within 0.50 Mpc of the clus- 
ter center and with m* — 1 + k s < m r < m* + 2 + k e . 
The fc-corrections are applied to individual galax- 
ies at a later stage, so these limits guarantee that 
we select all galaxies that can fall within m* — 1 < 
m r < m* + 2. The number of galaxies selected in 
the cluster region is N c ;„. 

2. We estimate the background contribution locally. 
We randomly select ten 20' x 20' boxes (avoiding 
bad areas) in a 1.3°-wide annulus, starting 3 h _1 
Mpc from the cluster center. Galaxies are selected 
within the same magnitude range as used for com- 
puting N c ; u . The median counts from the ten boxes 
are scaled to the cluster area to generate the back- 
ground estimate (Nftfc ff ). We adopt the interquartile 
range (IQR, which is the range between the first 
and third quartiles) as a measure of the error in 
Nfcfc s , which we term Qa-,bkg- The background cor- 
rected cluster counts (N c ;„ - Nf,^) is called N corr . 

3. Next, a bootstrap procedure is used to statistically 
apply k— corrections to the galaxy populations in 
each cluster. In each of 100 iterations, we ran- 
domly select N corr galaxies from those falling in the 
cluster region (N c /„). An elliptical k— correction 
is applied to 80% of the N corr galaxies, while an 
Sbc k— correction is applied to the remaining 20%. 
Finally, we use these k— corrected magnitudes to 
count the number of galaxies with m* — 1 < m r < 
m* + 2. The final richness estimate N ga ; s is given 
by the median counts from the 100 iterations. The 
richness error from the bootstrap procedure alone is 
given by Q a ,boot- The richness error is the combina- 
tion of this error and the background contribution, 

so that Q a = ^Ql Mot + Ql Mg - At this point, 

the richness error includes contributions from the 
fc-correction and the background galaxy correction, 
but not the redshift uncertainty, which is incorpo- 
rated in Step 5. 

4. If the cluster is too nearby or too distant, either 
the bright (m* - 1 + ks) or faint (m* + 2 + ke) 
magnitude limit, respectively, will exceed one of 
the survey limits (15.0 < m r < 19.5). We then 
apply the appropriate incompleteness correction to 
the richness estimate: 

r"V+ 2 $( m W m 
J 15 9(m)dm 



L* r -i *("*)<im 

72 = 720 — zr~rr~ (3) 

J „, <£>(m)dm 

We call 71 and 72 the low and high magnitude limit 
correction factors. 

5. The above steps are repeated using each of the 
n SUC cessfui photometric rcdshifts. The final clus- 
ter richness is recorded as the mean of the corre- 
sponding n SU ccessfui richnesses. We compute the 
mean of the richness errors from Step 3, as well 
as the dispersion among the richnesses from the 
^successful iterations. The former quantifies how 
much richness variation we expect based solely on 
cosmic variance, assuming there is no error in the 
redshift estimates, while the latter reflects the rich- 
ness error due to scatter in the photometric red- 
shifts. Because these are two independent sources 
of error, we add them in quadrature to derive the 
final richness error Q a ,Ngais- 

Total r-band luminosities (in solar units) and their er- 
rors are computed similarly to the richnesses. No at- 
tempt is made to fit and integrate luminosity functions 
for the individual clusters. 

2.5. Completeness and Contamination 

The globa l cont amination rate is estimated following 
lLopes et al.l ((2004). For each plate, we use the Raleigh- 
Levy (RL) distribution to generate N rea i t i x,y coordi- 
nates, where N rea i^ is the number of galaxies in the 
DPOSS catalog for plate i. Each galaxy in the RL cata- 
log is assigned a magnitude selected randomly from the 
real data. The density mapping and cluster detection 
is performed on the RL catalog for each plate, and the 
detected 'clusters' are assigned photometric redshifts at 
random from the real clusters in that plate. To esti- 
mate the global contamination rate for our sample, we 
ran this procedure on a plate-by-plate basis for the entire 
NoSOCS area. The results are shown in Figured! where 
the top panel shows the distribution of real clusters (solid 
line) and false clusters (from the RL simulations, dotted 
line), as a function of richness. The total contamination 
rate is 8.4%, consistent with our goal of achieving < 10% 
contamination when setting the cluster detection param- 
eters. The bottom panel shows the contamination rate as 
a function of richness. For very rich clusters (N ga i s > 50) 
the contamination rate is negligible, and only rises above 
5% for N gaU < 20. 

The redshift- and richness-dependent completeness 
functions for each plate are provided in Table [T] The 
first column gives the plate number; for each plate, 
there are 42 entries, using six richnesses (N ga i s = 
15,25,35,50,80,120) given in the second column, at 
seven redshifts (z = 0.08 to 0.32 with 5z = 0.04) given 
in the third column. The fourth column gives the recov- 
ery rate (in percent) of clusters with the given richness, 
at the listed redshift, for that specific plate. To use a 
large cluster sample for cosmology requires knowledge of 
the mass-dependent selection function (SF). Currently, 
there are two methods to generate these, either by em- 
pirically calibrating a mass-observable relation, or using 
large simulations to construct mock galaxy catalogs from 
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TABLE 1 

Completeness Functions 



Plate 


Richness 


Redshift 


Completeness 


005 


015 


0.08 


34.0 


005 


015 


0.12 


10.0 


005 


015 


0.16 


10.0 


005 


015 


0.20 


2.0 


005 


015 


0.24 


0.0 


005 


015 


0.28 





005 


015 


0.32 


0.0 


005 


025 


0.08 


78.0 


005 


025 


0.12 


68.0 


005 


025 


0.16 


50.0 


005 


025 


0.20 


12.0 


005 


025 


0.24 


8.0 


005 


025 


0.28 





005 


025 


0.32 





005 


035 


0.08 


96.0 


005 


035 


0.12 


82.0 


005 


035 


0.16 


72.0 


005 


035 


0.20 


30.0 


005 


035 


0.24 


8.0 


005 


035 


0.28 


8.0 


005 


035 


0.32 





005 


055 


0.08 


98.0 


005 


055 


0.12 


92.0 


005 


055 


0.16 


92.0 


005 


055 


0.20 


82.0 


005 


055 


0.24 


46.0 


005 


055 


0.28 


10.0 


005 


055 


0.32 


6.0 


005 


080 


0.08 


100.0 


005 


080 


0.12 


98.0 


005 


080 


0.16 


94.0 


005 


080 


0.20 


92.0 


005 


080 


0.24 


64.0 


005 


080 


0.28 


34.0 


005 


080 


0.32 


8.0 


005 


120 


0.08 


100.0 


005 


120 


0.12 


100.0 


005 


120 


0.16 


98.0 


005 


120 


0.20 


100.0 


005 


120 


0.24 


90.0 


005 


120 


0.28 


64.0 


005 


120 


0.32 


24.0 



which clusters are selected. iKoester et al.l (|2007al) use 
the latter to estimate the purity and completeness of the 
SDSS MaxBCG cluster catalog, but we do not have such 
simulations corresponding to our data. Using X-ray ob- 
servations, we are in the process of developing an op- 
timized richness estimator to generate an N ga [ s — Mass 
relation. However, that work is beyond the scope of this 
paper. 

2.6. Other Changes 

Following Paper II, we generate ten additional AK 
maps for each plate, using a set of galaxy catalogs 
for each plate with random photometric zero-point off- 
sets added to the r-band magnitudes, drawn from the 
k nown photometri c error distribution for DPOSS given 
in iGal et a l. ( 2004) . In Paper II we required that a clus- 
ter candidate be detected in seven of the ten zero-point- 
error-added maps in addition to the original map. In the 
final catalog, we now require only that a candidate be 
detected in any seven of the eleven maps, as there is no a 
priori reason to give preference to the original map. This 
typically results in ~ 5 — 10% additional candidates per 
plate. 

3. CLUSTER MORPHOLOGICAL PROPERTIES 
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Fig. 4. — Contamination of the entire NoSOCS catalog by false 
clusters. The top panel shows the distribution of real clusters (solid 
line) and false clusters (from the RL simulations, dotted line), as a 
function of richness. The total contamination rate is 8.4%, consis- 
tent with our goal of achieving < 10% contamination when setting 
the cluster detection parameters. The bottom panel shows the 
contamination rate as a function of richness. For very rich clusters 
(N ga i s > 50) the contamination rate is negligible, and only rises 
above 5% for N gats < 20. 

3.1. Substructure Measures 

Four substructu re measures are computed for each 
candidate cluster (jLopes et al.l 12006). Only clusters at 
0.069 < z p hot < 0.196 are examined; in this redshift 
range, we completely sample the cluster luminosity func- 
tion spanning m* — 1 < m r < m* + 1. We apply the 
angular separation test (AST), the Fourier elongation 
test (FE), the Lee statistic (Lee 2D), and the symme- 
try test (/3) to 10,575 clusters, within a radius of 1.5/i _1 
Mpc around the recentered positions, and a significance 
level threshold of 5%. The rationale for these choices are 
discussed in §5 of lLopes et a l. (2006), wh ile detailed de- 
script ions of all four tests are provided bv lPinknev et al.1 
( 1996). Very briefly, the values taken on by the four tests 
indicate substructure as follow: 

• (3: For a symmetric distribution < (3 >« 0, while 
values of <(3> greater than indicate asymmetries. 

• AST: This statistic takes on values near unity for 
substructure- free systems, and less than 1.0 for 
clumpy distributions. 

• FE: Values of this statistic greater than 2.5 indicate 
significant deviations from circularity. 

• Lee 2D: Larger values of this statistic indicate the 
presence of two subclumps in the galaxy distribu- 
tion. 
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TABLE 2 

Statistics From Three NoSOCS Regions 



Region N c!uaters Area {deg 2 ) p (N deg 2 ) z med z QtT N galStmed N galStQa P(z phot ) P(N gals ) 

NGP, good 7985 5681.31 1.405 0.1416 0.057 18.65 9.73 

NGP, poor 3491 2812.54 1.241 0.1442 0.057 18.96 9.64 0.259 0.326 

SGP 4026 2917.28 1.380 0.1309 0.061 17.83 9.72 0.000 0.312 

Combined 15502 11411.13 1.358 0.1393 0.058 18.50 9.75 



TABLE 3 
Cluster Properties 



Name 


RA 


Dec 


Ndet 


< Zphot > 




Ngals 




Lopt 


v{L op t) 


<off> 


& 


NSC J000016+103643 


0.06709 


10.61203 


11 


0.1319 


0.0109 


18.2 


4.9 


0.468 


0.173 


0.17 


29.50 


NSC J000018+204800 


0.07902 


20.80006 


10 


0.0901 


0.0060 


8.7 


4.2 


0.113 


0.102 


0.10 


-33.80 


NSC J000020+210327 


0.08433 


21.05770 


7 


0.1674 


0.0014 


23.0 


3.6 


0.647 


0.230 


0.16 


-21.80 


NSC J000024+ 142904 


0.10351 


14.48461 


11 


0.1254 


0.0051 


25.5 


5.7 


0.574 


0.167 


0.10 


50.20 


NSC J000029+215512 


0.12189 


21.92005 


11 


0.1332 


0.0035 


30.1 


1.0 


0.550 


0.152 


0.05 


-13.10 


NSC J000032+141432 


0.13691 


14.24238 


11 


0.0839 


0.0150 


12.1 


6.1 


0.302 


0.152 


0.48 


-23.50 


NSC J000038+063046 


0.16178 


6.51303 


11 


0.2276 


0.0063 


33.8 


4.3 


1.201 


0.565 


0.42 




NSC J000040+065659 


0.16740 


6.94999 


7 


0.1389 


0.0038 


19.6 


4.0 


0.383 


0.149 


0.09 


38.40 


NSC J000048+125623 


0.20284 


12.94000 


11 


0.1087 


0.0044 


15.1 


4.1 


0.234 


0.109 


0.27 


53.10 


NSC J000051+152013 


0.21461 


15.33697 


11 


0.1254 


0.0044 


18.6 


6.3 


0.318 


0.139 


0.35 


66.90 


NSC J000056+004551 


0.23454 


0.76427 


11 


0.2114 


0.0298 


34.5 


3.5 


1.186 


0.377 


0.27 




NSC J000057+064615 


0.24058 


6.77092 


7 


0.2118 


0.0022 


0.0 


0.0 


0.000 


0.000 


0.09 




NSC J000105+023236 


0.27109 


2.54333 


11 


0.0927 


0.0017 


18.5 


4.2 


0.420 


0.147 


0.13 


16.00 


NSC J000125+181149 


0.35463 


18.19711 


11 


0.1185 


0.0020 


16.4 


3.6 


0.262 


0.122 


0.10 


-19.30 



TABLE 4 

Substructure Measurements 





Name 


P 


FE 


LEE2D 


AST 


NSC 


J062054+861617 


-14.9 


2.243 


2.667 


33.3 


NSC 


J074613+854032 


-43.2 


1.673 


1.644 


21.0 


NSC 


J054111+843927 


-15.1 


1.413 


1.547 


25.6 


NSC 


J065439+845907 


-29.9 


2.798 


1.722 


17.2 


NSC 


J054209+842633 


18.0 


0.989 


1.623 


12.8 


NSC 


J055822+841733 


-38.6 


1.560 


2.332 


17.8 


NSC 


J061210+841036 


-35.4 


1.629 


1.650 


32.7 


NSC 


J065407+842104 


-32.4 


0.485 


1.968 


19.5 


NSC 


J064833+841519 


84.5 


1.180 


2.259 


14.8 


NSC 


J073249+841701 


1.2 


0.913 


1.676 


22.1 


NSC 


J093356+845601 


-17.5 


1.395 


1.296 


15.8 


NSC 


J094741+844440 


7.0 


0.991 


1.386 


21.8 


NSC 


J094540+843709 


2.8 


2.423 


2.214 


17.6 


NSC 


J083925+835412 


-15.5 


1.842 


2.180 


27.8 


NSC 


J083042+824948 


24.4 


2.217 


1.842 


23.1 


NSC 


J085200+830113 


-86.9 


0.425 


2.119 


27.5 


NSC 


J091312+825157 


17.8 


1.785 


2.086 


14.0 


NSC 


J084436+861547 


-1.9 


1.703 


2.389 


17.9 


NSC 


J105955+853131 


17.1 


2.475 


2.048 


31.2 


NSC 


J130353+844618 


-4.5 


1.241 


1.572 


26.4 


NSC 


J104433+840151 


111.6 


1.027 


13.896 


8.6 



The main data table (Table [3|) includes only the /3-test 
results, while Tab le [H provides the results of all four 
tests. As noted in lLopes et al.l (|2006l ). the /? test is the 
most sensitive to substructure. 

3.2. Estimating the Dynamical Radii 

We attempt to estimate the typical length scale charac- 
terizing the virializcd regions of the clusters of our sam- 
ple. Both the the ory of gravitational collapse in an ex- 
panding Universe IGunn k Gottl [19721 e.g. and N-body 
simulations suggest that the virialized mass of a clus- 
ter is generally contained inside the surface where the 
mean interior density is about 200 times the critical den- 



sity p c (z), at the redshift of the cluster (|Carlberg et al.1 
1997): 



(Pm), 



200p c (z) 



200 



A I 



Pm( z ) 



(4) 



where (pm)r is the mean mass density of the cluster 
within R and ~p M (z) is the mean mass density of the 
Universe at redshift z. We assume that the radial distri- 
bution of galaxies within a cluster follows the dark matter 
and neglect possible variations of the mean mass of galax- 
ies, rrigai, with environment. With these simplifications, 
(Spm/Pm) = (Spgai/Pgai) — (SVgai/Vgai), where we de- 
fine the number density of galaxies, v ga i = Pgai/rn-gai- 
Thus, from Eq. 4 we get a simple formula relating the 
number density to i?20o: 



/ \ 200 - ( \ 



A I 



(5) 



Finally, the spatial mean number density of galaxies 
Wgai) appearing in this formula may be related to the 
observed projected number density (S ga ;), through the 
approximation: {v ga i) R ~ vri? 2 (T, ga i} R /(47r/3)i? 3 = 

| (£ ga i} /R. 

We calculate -R200 within the same redshift range used 
for substructure measurements (0.069 < z < 0.196), but 
further limited to N ga i s > 25 (2681 clusters). We then 
select only those clusters with less than 10% of their 
area within a circle of radius 1.5/i _1 Mpc intersected 
by projected circles from neighboring clusters, leaving 
1637 clusters. The 10% overlap criterion avoids struc- 
tures whose projected profiles and background regions 
are likely contaminated by galaxies from a neighboring 
cluster. Applying our methods iteratively could be used 
to relax this criterion but we have not done so as the 
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Fig. 5. — -R200 measurements plotted against richness for the 
"best" sample. The data were binned, subject to a minimum num- 
ber of objects/bin, and medians were taken. The errors bars give 
the 1<t range of values for i?200 in each bin. The best fit relation 
between R200 and richness is shown. 

uncertainties on -R200 al ready significant. T he lu- 
minosity function derived bv lBlanton et a l. (2001), inte- 
grated to the completeness limit of the NoSOCS catalog, 
M r = —19.8, gives a good estimate of the mean num- 
ber density of galaxies in the Universe at redshifts < 
0.2. Since the NoSOCS counts are complete only in a 
restricted apparent magnitude range of 15.0 < r < 19.5, 
for each cluster we computed a completeness correction 
factor considering the absolute magnitude limit above. 
T hese were estimated u sing the luminosity function given 
by iPaolillo et al.l (|200lD . which was derived from a sam- 
ple of Abell clusters detected in the DPOSS survey. 

Foll owing the procedure described by lLopes et al.l 
(2006), for each cluster, the background density contri- 
bution was calculated using an annular ring about the 
cluster center with inner and outer radii of Ri„ — 3 Mpc 
and R ou t = 4.6 Mpc, respectively. The cumulative pro- 
jected number density profiles appearing in Eq. 5 are 
then calculated by counting galaxies in concentric annuli 
around the cluster center. The ring widths are variable, 
defined by requiring a constant number of galaxies per 
ring. These counts were then corrected for the back- 
ground contribution and for completeness. Further cor- 
rections were applied to account for the regions within 
the annuli that intersected the bad areas due to bright 
objects or the densitometry spots. Furthermore, when 
computing i?200 for a given cluster, areas around neigh- 
boring clusters were masked with a 1.5 Mpc radius to 
avoid projection effects, resulting in large excluded ar- 
eas for low redshift clusters. For each cluster, its galaxy 
number density, was calculated by excluding galaxies lo- 
cated in the overlap areas and correspondingly correcting 
the counting areas. For each cluster, the ratio f overlap of 
the number of galaxies in the overlapping areas to the 
total number of galaxies within a maximum search circle 
centered on the cluster center (R se arch), was estimated. 
Clusters with f overlap > 0.7 did not have -R200 computed. 
For 17 of 1637 clusters the measurement failed because 
the computed values were unphysically large (-R200 > 4 



Mpc h^ , thus extending to the background area). For 
2 clusters the background density was too high and no 
meaningful density profile could be obtained. 

The solution i?2oo of Eq. 5 is obtained by spline in- 
terpolating the cumulative density profile using the 5 
points nearest to the solution. Figure [5] show results 
for the "best" clusters as a function of richness. This 
subsample consists only of clusters with high richness 
(Ng a i s > 25), chosen to reduce the effects of background 
fluctuations. Furthermore, clusters whose analysis re- 
gions were affected by bad areas or neighboring clusters 
over > 50% of their total projected areas were discarded, 
as were those which crossed plate boundaries. 

Examination of the -R200 values shows that they spa n 
the same range as those reported in lHansen et al.l |2005), 
but we find very large scatter as a function of N ga i . This 
is seen in the error bars in Fig. which are large de- 
spite having ~ 20 clusters per bin. This is likely due 
to a combination of shallow depth, large photometric er- 
rors and the exclusion of significant regions due to cluster 
overlaps. Nevertheless, the overall relation between -R200 
and Ng a i s is reasonable, as shown in Fig.O A linear best 
fit to this data with R 200 oc N£ al yields a = 0.40 ± 0.02, 
well wit hin the expe c tation s from the results of the anal- 
ysis by lLopes et al.l (|2006| ) . In that work it was shown 
that for X-ray clusters in common with a subsample of 
NoSOCS clusters without substructure, Tx oc N^ al , with 
~ 0.8. Since M(R 20 o) oc T 3 / 2 , where M(R 200 ) is the 

3 /2 

cluster mass inside -R200, it follows that i?200 oc , as 
we have found here. 

These findings are comparable to those in the litera- 
ture. The range of i?2nn spanned by our clusters is sim- 
ilar to those in lHansen et alj (|2005f ) when transformed 
to their cosmology with h = 1, although we do not ex- 
tend to the lowest richness systems (N ga i Si MaxBCG < 4) 
that they include. They also find, using a richness mea- 
sured solely from the red sequence in the MaxBCG tech- 
nique, a = 0.57 in the r' band. Similarly. iPopesso et alJ 
(|2007f l. examining clusters detected in both the RASS 
and SDSS, find N 2 oo oc M^ 1 ; ass uming mass scales 
with volume this yields a = 0.37. ICollister fc Lahavl 
(2005) looked at clusters and groups in the 2dFGRS 
and derived N oc M°- 99 , whic h gives a — 0.34. Using 
K-band data, iLin et all (|2004l ) find N gal cx M^f, or 
a = 0.39. The aforementioned surveys compare rich- 
ness and mass both measured within i?200- A direct 
comparison -R200 to richne ss measured in a fixed phys- 
ical aperture was done bv lYee fe Ellingso"n1 d2003j), who 
used a radius of 0.5 Mpc /ijjj, 1 when examining CNOC 
clusters; they find -R200 oc B^f 7 . As noted above, we 
find -R200 oc ^gais' brooby consistent with all of these 
results despite the large differences in richness measure- 
ment techniques. The photometric data, cluster detec- 
tion and especially richness measurements are all dis- 
tinct, and a full comparison would require running analo- 
gous detection and richness codes on both datasets. Such 
bidirectional tests will be fundamental to assessing sys- 
tematic effects in cluster catalogs. 

4. GLOBAL SAMPLE PROPERTIES 
4.1. Comparison of Three Regions 
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- NGP well calibrated 

NGP poorly calibrated 
(scaled) 
SGP (scaled) 
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- NGP well calibrated 

NGP poorly calibrated 

(scaled) 

SGP (scaled) 



Fig. 6. — Redshift distributions for the three independent regions 
of NoSOCS.. The NGP poorly calibrated (dotted line) and SGP 
(dashed line) distributions are scaled to the same total number of 
clusters as in the well-calibrated NGP region (solid line). 



Fig. 7. — Richness distributions for the three independent regions 
of NoSOCS. The NGP poorly calibrated (dotted line) and SGP 
(dashed line) distributions are scaled to the same total number of 
clusters as in the well-calibrated NGP region (solid line). 



As discussed earlier, the cluster catalogs presented here 
and in Paper II are generated for three independent re- 
gions. Although the photometric calibration, object clas- 
sification, and cluster detection are all performed in an 
identical manner, one may still expect systematic vari- 
ations between these areas, especially when considering 
the poorly calibrated NGP region. With the extremely 
large number of cluster candidates in the three regions, 
we expect that the redshift and richness distributions 
should be very similar. We show the results of this com- 
parison for redshifts in Figure [6] and richnesses in Fig- 
ure [7] The histograms for all areas are scaled to the same 
total number of clusters as the well-calibrated NGP re- 
gion. Table [5] gives the region name, number of clusters, 
total area and projected density of clusters in Columns 
1 through 4, as well as the median and Q a for z p \- lot 
(Columns 5 and 6) and N ga i s (Columns 7 and 8). The 
P- values from Kolmogorov-Smirnov tests comparing the 
redshift and richness distributions for both the poorly 
calibrated NGP region and the SGP region to the well- 
calibrated NGP region are provided in Columns 9 and 
10 of Table [5] We test the redshift distributions over 
the range 0.07 < z p hot < 0.3 (where the completeness is 
high) and find that they are consistent for the two NGP 
areas, while the SGP region is discrepant, with an ex- 
cess of clusters at z < 0.13. Beyond this redshift, the 
SGP redshift distribution agrees very well with the other 
regions. We also examined the contamination rates esti- 
mated in §2.4 separately for each region. The results are 
shown in Figure [8] The three regions are evidently very 
similar, although the poorly calibrated NGP area may 
be slightly worse, as we would expect with less accurate 
photometry. Whether the differences can be attributed 
to cosmic variance or not is unclear; structures on scales 
of ~ 300/i _1 Mpc, corresponding to a redshift range lo- 
cally of Az ~ 0.1 are seen in 2MASS a nd other surveys 
(jFrith et alj|2003t ICourtois et al.ll200l . 




Fig. 8. — Contamination rates in the three separate NoSOCS re- 
gions. From left to right, we show the well-calibrated NGP, poorly 
calibrated NGP, and SGP, while the top and bottom rows are the 
same as Fig. [4] 

The complete catalog of 15,502 clusters is presented in 
Table [31 The columns in this table are: 

1. Cluster Name. The name is NSC (for Northern 
Sky Cluster), followed by the coordinates JHH- 
MMSS+DDMMSS. 

2. Right ascension in J2000.0 decimal degrees. For 
clusters where the photometric redshift estimator 
succeeded, this is the mean of the recentered posi- 
tions. Where the photo- z failed, this is the original 
detected position. 

3. Declination in J2000.0 decimal degrees. See notes 
for RA. 

4. The number of times ridet this cluster was detected 
in the 11 detection passes (see §2.5). 



4.2. The Final Catalog 



5. The mean photometric redshift, z p hot, from the ten 
photo- z runs. 
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6. The photometric redshift error, including the con- 
tribution from the scatter in the photo-z relation 
and the multiple photo-z runs. 

7. The mean richness N ga i s from the ten richness runs. 

8. The richness error, including contributions from 
the ^-corrections, background variance, and red- 
shift errors. 

9. The r-band optical luminosity L op t, in solar units. 

10. The luminosity error. 

11. The (3 substructure parameter. This was only cal- 
culated for clusters at 0.069 < z p hot < 0.196. 

12. The mean offset (in Mpc) from the original de- 
tected position in the ten photo-z runs. If the 
photo-z failed, this is left blank. 

For the subset of 2,681 clusters with N ga i s > 25 and 
0.069 < z < 0.196, we provide the X-ray luminosities 
measured within fixed apertures of 0.5 and 1.0 hr x Mpc 
(in units of 10 43 erg s -1 ) along with the associated errors 
and X-ray temperatures in Table The derivations of 
the X-ray quantities are discussed in §6 below. Table 
131 provides the results of all four substructure tests for 
10575 clusters at 0.069 < z < 0.196. 

5. COMPARISON TO THE SDSS MAXBCG CATALOG 

It is instructive to compare large cluster catalogs cover- 
ing the same sky area, both as a consistency check for the 
newer catalog and to search for possible systematic er- 
rors. In our earlier work we compared the first NoSOCS 
area to the Abell catalog. Since then, a new, deeper 
cluster catalog based on Sloan Digital Sky Survey data 
and using the MaxBC G algorithm has been published by 
iKoester et al.l (|2007bh . In this section we compare our 
catalog to theirs, examining recovery rates and richness 
estimates. 

The Sloan Digital Sky Survey (jYork et alj I2000D has 
been used to generate a variety of cluster catalogs with 
different techniques , some of which are compared in 
iBahcall et al.l (|2003f ). However, that work used only a 
small area of the sky with early SDSS, spanning a few 
hundred square degrees. The only cluster detection tech- 
nique applied to a majority of the SDSS sky coverage 
is the red sequence - brightest cluster gal axy technique 
called MaxBCG, described in de tail in IKoester et al.l 
(|2007aft . The sample described in (jKoester et al.ll2007bft 
covers ~ 7500 deg 2 , containing nearly 14,000 clusters. 
The increased depth, higher photometric accuracy and 
multiple passbands of SDSS allow for the generation of 
cluster catalogs that are more complete for poor clusters, 
extend to higher redshifts, and yi eld better photometric 
redshift estimates. Furthermore, IKoester et al.l (|2007af ) 
have used cosmological simulations to assess their com- 
pleteness and false detection rate as a function of cluster 
mass. Thus, their work can be used as a benchmark for 
the NoSOCS sample, which covers a larger sky area. 

It is important to remember that the bright flux limit 
of the galaxy catalog used here makes NoSOCS an es- 
sentially flux-limited sample. This can be seen in Fig. 
3 of Paper II; the completeness of our survey is highly 
richness dependent even at z ~ 0.2. In contrast, the 



SDSS photometric catalog is ~ 3 magnitudes deeper. 
The MaxBCG method relies on the E/SO ridgeline to 
detect clusters, and samples such galaxies down to 0.4X» 
out to z = 0.4. Thus, the MaxBCG catalog, trimmed 
to z = 0.3 to reduce photometric redshift uncertain- 
ties, provides something close to a volume-limited sam- 
ple. The completeness is near unity for all cluster masses 
>3_x 1O 14 M , as shown in Fig. 7 of IKoester et all 
(2007tJ). Only for poor systems is this untrue, since the 
limit of N ga - St MaxBCG > 10 imposed on the published 
sample will introduce some incompleteness. The distinc- 
tion between our flux-limited catalog and the volume- 
limited MaxBCG catalog is evident in the mutual recov- 
ery rates discussed below. 

First, we checked the SDSS sample for clusters falling 
into any of our bad areas. Only 141 of their 13,823 
clusters (1.02%) are eliminated in this way. This sug- 
gests that the sizes of our exclusion regions are reason- 
able, especially since the long exposures on the pho- 
tographic plates yield larger saturated regions around 
bright stars than the shorter SDSS exposures. We re- 
strict our comparison to the NGP region bounded by 
135° < RA < 225°, 0° < Dec < 50° since the SDSS data 
covers only small strips in the SGP and does not extend 
to the northernmost declinations. A rectangular region 
is also trimmed from both catalogs to account for a miss- 
ing stripe in the SDSS area. These cuts result in an over- 
lap region of 3100 deg 2 containing 5,595 maxBCG and 
4,275 NoSOCS clusters. Further restricting our catalog 
to the same redshift range as MaxBCG (0.1 < z < 0.3) 
leaves only 3,299 NoSOCS clusters. However, applying 
these cuts based on our noisy photometric redshift esti- 
mator will in troduce compl e x effe cts in the comparisons, 
as noted by IBahcall et al.l (|2003l ). so we do not apply 
this cut to our catalog. Figure [H] shows the region of 
sky used, with NoSOCS clusters in the top panel and 
maxBCG clusters in the bottom panel. While a gener- 
ally good correspondence is seen, there are clearly many 
clusters not in common to the two catalogs. The overall 
large scale structure, including filaments spanning tens 
of degrees, is well reproduced by both surveys. 

5.1. Does MaxBCG Find NoSOCS Clusters? 

We first examine the recovery rate of our clusters in 
the MaxBCG catalog. We simply search for the clos- 
est projected match to each NoSOCS cluster among the 
SDSS clusters. The angular separation between matched 
clusters is converted to a physical distance in kilopar- 
secs using the NoSOCS photometric redshift. Of the 
4,275 NoSOCS clusters, only 49.3% are matched to a 
MaxBCG counterpart within 1 Mpc. The top panels of 
Figure [10] show the recovery rate of NoSOCS clusters by 
MaxBCG as a function of matching radius and NoSOCS 
richness. For poor clusters (N ga i Sy NoSOCS < 30) the 
recovery rate is low, even using large matching radii. 
This suggests that (a) MaxBCG may fair poorly at de- 
tecting poor systems which have weak or no red se- 
quence and no BCG, and/or (b) the contamination rate 
in our catalog is high for poor clusters. On one hand, 
the MaxBCG catalog demonstra tes a completeness o f 
> 80% for N galsM axBCG > 10 (|Koester et al .l l2007alh 
based on both Monte Carlo simulat ions where Abell-type 
clusters are inserted into the data (jKoester et al.l l2007b) 
as well as cluster detection run on large mock catalogs 
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Fig. 9. — Comparison of the projected distribution of clusters in our catalog (top) and the SDSS MaxBCG catalog (bottom), 
excluded rectangular area corresponds to a stripe missing from the MaxBCG catalog. 



The 



(jRozo et al.ll2007j ). It is also nearly volume- limited and 
should therefore contain all such structures at z < 0.3; 
however, the a posteriori limit of N ga 



gals, MaxBCG 



> 10 



imposed on the published catalog is likely to have elimi- 
nated many poor systems that we detect. On the other 
hand, based on Figure [4] only ~ 15% of poor NoSOCS 
clusters are expected to be false detections. The low re- 
covery rate of poor NoSOCS clusters by MaxBCG calls 
into question either one or both of these results, and 
requires further detailed study, especially using spectro- 
scopic redshifts to determine the reality of these systems. 
We suspect that many of the poor systems we detect 
but not in the published MaxBCG catalog may have 
MaxBCG richnesses below their publication threshold. 

More intriguing is the low ~ 75% recovery rate for 
very rich NoSOCS clusters. Our estimated co mpleteness 
is ~ 95% for N gals > 50 out to z = 0.25 (|Gal et all 
l2003h . while the con tamination rate is neg ligible for rich 
clusters. Similarly, iKoester et al.l (|2007af) claim nearly 
100% completeness for similar systems. Examination of 
the unrecovered systems shows that they are typically 



at 0.2 < Zphot < 0.3, as shown in the top panel of Fig- 
ure [TlJ This suggests that a combination of the strict a 
posteriori redshift limits imposed on the MaxBCG cat- 
alog, along with the significant scatter in the NoSOCS 
photometric redshifts is responsible for a significant por- 
tion of the observed incompleteness. It is unlikely that 
richness errors cause this incompleteness, since the two 
surveys' richnesses are well correlated. Nevertheless, it 
will be important to carefully examine the rich systems 
found by only one of the techniques to understand po- 
tential biases. These comparison difficulties also show 
that applying a posteriori limits (in richness, redshift, or 
some other property) to publicly available catalogs makes 
them troublesome (if not impossible) to use in such com- 
parative studies. 

5.2. Does NoSOCS Find MaxBCG Clusters? 

Next, we reverse the sense of the comparison, exam- 
ining the completenes s of our catalog relative to that of 
IKoester et~aTJ (|2007bh . Here, we use the MaxBCG cat- 
alog as the fiducial source, searching for the nearest (in 
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Fig. 10. — Top: Reco very rate of N oSOCS clusters in the 
MaxBCG catalog of Koe ster et al.1 l(2007b) as a function of match- 
ing radius. Bottom: The reverse comparison, showing the recov- 
ery of MaxBCG clusters in the NoSOCS catalog. 
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Fig. 11.— Top: The redshift distribution of NoSOCS clusters 
not found by MaxBCG. The different line types denote distinct 
intervals of NoSOCS richness. Bottom: The redshift distribution 
of MaxBCG clusters not identified in NoSOCS. The different line 
types denote distinct intervals of MaxBCG richness. 



projection) NoSOCS cluster. Angular separations are 
converted to physical distances using the MaxBCG pho- 
tometric redshifts. Because the MaxBCG catalog should 
be essentially 100% complete in the redshift range probed 
by NoSOCS, it provides a potential basis for testing our 
own completeness (but see the caveats above). The re- 
sults are shown in the top panel of Figure [TU] with the 
recovery rate of MaxBCG clusters by NoSOCS as a func- 
tion of matching radius and MaxBCG richness. It is im- 
mediately apparent that NoSOCS does extremely well at 
discovering rich clusters, finding 80-100% of the richest 
MaxBCG clusters. Measured this way, NoSOCS is more 
complete for the richest clusters than MaxBCG, although 
this may be due to clusters falling below the z > 0.1 limit 
imposed on the MaxBCG catalog. However, for poor 
clusters and groups, the recovery rate is low. Nearly half 
of the MaxBCG clusters have N ga i St MaxBCG < 15, falling 
into the lowest richness bin in this plot. Clearly, neither 
our algorithm nor that of MaxBCG is anywhere near 
complete for group-mass systems. The bottom panel of 
Figure [11] shows that the NoSOCS completeness drops 
with redshift to z ~ 0.3, as expected for our flux-limited 
survey. Nevertheless, the completeness is very high even 
for moderately poor systems to z ~ 0. 2, consistent wit h 
the estimates shown in Fig. 5 and 6 of lGal et alj |2003). 

5.3. Comparison of Cluster Properties 
5.3.1. Photometric Redshifts 

Beyond examining the completeness of these cluster 
catalogs, we use the more accurate MaxBCG photomet- 
ric redshifts to test our own estimates. We also exam- 
ine the relationship between the NoSOCS and MaxBCG 
richnesses. The left panel of Figure [12] shows the com- 
parison of photometric redshift estimators, as a func- 
tion of NoSOCS richness, for NoSOCS clusters with a 
MaxBCG counterpart within 0.75 Mpc. The poorest 
clusters are shown as the smallest dots. Open squares 
show the MaxBCG photometric redshifts on the ordi- 
nate, and their spectroscopic redshifts on the abscissa. 
For poor clusters (small dots), the scatter between the 
two estimators is high, and the NoSOCS z p hot appears 
to underestimate the true redshift. For clusters with 
N ga is 7 NoSOCS > 20, the scatter is dramatically reduced 
and there is only a small offset, which disappears for 
Ngais, NoSOCS > 30 . These are quantified in the right 
panels of Figure [TSJ which shows the scatter (top) and 
median offset (bottom) between the two photometric red- 
shift estimators, as a function of N ga i St N socs an d the 
matching radius. Assuming the MaxBCG measurements 
are more accurate, we overestimate the redshifts of poor 
clusters. This may be due to the training sample used, 
which consists almost exclusively of Abell clusters, which 
are much richer than these poor groups. Furthermore, 
because of the minimum redshift (z > 0.1) imposed on 
the MaxBCG sample, there is a bias at the low redshift 
end, where most of the poor NoSOCS clusters are de- 
tected. In fact, there are only two clusters in the sample 
with N gals , No socs < 10 and 0.15 < z phot:N oSOCS < 0.25 
matched within 750kpc. If we move to the next richness 
bin, 10 < Ng a i Si NoSOCS < 20, there are 115 clusters 
with 0.15 < Zphot, NoSOCS < 0.25, of which 80% are at 
z < 0.19. This effect is shown by the line and aster- 
isk in the bottom right panel of Figure [T2l where the 
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Fig. 12.— Left: Comparison of NoSOCS and MaxBCG rcdshifts 
for clusters matched within 0.75Mpc. The solid circles show the 
NoSOCS photometric redshifts compared to those from MaxBCG. 
The point size increases with increasing richness. Open squares 
show the MaxBCG spectroscopic redshifts vs. their photometric 
redshifts. Right: The dispersion (top) and median offset (bot- 
tom) between NoSOCS and MaxBCG photometric redshifts, as a 
function of richness and matching radius. 

median and scatter of the z p h ot differences is computed 
only for those clusters with 10 < N ga i B ^ No socs < 20 and 
0.15 < z p hot, NoSOCS < 0.25, compared to the solid circle 
if no rcdshift cut is used. Applying this limited redshift 
range reduced the median offset by ~ 50%. 

5.3.2. Richness 

The number of galaxies in a cluster may be directly 
related to the underlying dark halo mass. If this is true, 
purely photometric cluster surveys are adequate to con- 
struct the cluster mass function, and in concert with pho- 
tometric redshifts, measure its evolution. 

To test the reliability of such richness estimates, we 
compare our richnesses to those from the MaxBCG cat- 
alog. Both surveys compute richnesses in fixed phys- 
ical apertures as well as within T2oo- The results are 
show in Figure [13j where the left panel shows our N ga i 8 
(within a 500kpc radius aperture) versus MaxBCG rich- 
nesses in the same aperture, and the right panel com- 
pares our fixed-aperture richness with MaxBCG's i?200 
richness. We only use clusters whose centroids agree to 
within 500 hr 1 kpc between the two surveys, and with 
0.1 < Zpf lo t < 0.3, resulting in a sample of 1,072 clus- 
ters. The small blue points show all matched clusters, 
while red points are those matches where the NoSOCS 
and MaxBCG photometric redshifts differ by less than 
0.03. The large open squares show the medians in bins 
of width ANg a i s — 10, along with the rms scatter. Al- 
though there is a moderately large dispersion between 
our richnesses and those from MaxBCG, they are well 
correlated. The solid lines in Figure [TBI show the best-fit 
relations: 



N, 



gals, Max 



bcg = 12-55 + 0.260 x N g 



(6) 



and 



N R200 ,MaxBCG = 8.60 + 0.375 X Ngai 



(7) 

As expected, the richnesses computed within the same 
fixed apertures (0.5 Mpc) are much better correlated. 
However, these relations should not be used to convert 
between richnesses from the two surveys. The MaxBCG 
catalog is censored at low richness (where most of the 
clusters are found) and the scatter in the relation is very 
high. Much of the scatter is likely due to the differ- 
ent definitions of richness, where we count all galaxies 
while MaxBCG effectively counts only galaxies along the 
E/S0 ridgeline. Due to the large photometric errors in 
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Fig. 13. — Comparison of NoSOCS and MaxBCG richnesses for 
clusters with 0.1 < z p hot < 0-3 matched within 0.5Mpc. Small blue 
points show all of these clusters, while larger red points require that 
the photometric redshifts agree within Az < 0.03. Large square 
show the medians in bins of width AN ga i a = 10. The left panel 
compares richnesses measured in the same fixed aperture of 500 
kpc. The right panel compares our fixed-aperture richness with 
the MaxBCG richness within their -R200 ■ The best fit relations are 
shown as the solid lines. 
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Fig. 14. — NoSOCS photometric redshifts vs. NORAS spec- 
troscopic redshifts. Points are color coded by the closeness of the 
positional match. 

DPOSS, we cannot replicate a richness using only the 
red-sequence galaxies. Additional scatter is introduced 
by the different photometric redshifts changing the an- 
gular sizes of the apertures between the two surveys. 

6. X-RAY MEASUREMENTS 

6.1. Comparison to NORAS 

It is instructive to compare the results of optical and 
X-ray cluster surveys, for the purposes of examining com- 
pleteness and testing properties (N ga is, Lx, L opt ) that 
might be useful as mass proxies. The largest existing 
X-ray survey in the Northern hemisphere is NORAS 
(|B6hringer et al.ll2000f ). with 378 clusters at 5 > 0°. NO- 
RAS is useful not only as the largest, homogeneous cata- 
log of low-to-moderate redshift clusters, but also because 
spectroscopic redshifts have been obtained for the entire 
sample. To match this sample to NoSOCS, we first re- 
move 55 NORAS clusters in our bad areas, and an addi- 
tional 97 at low galactic latitude. This leaves a sample 
of 226 NORAS clusters for comparison, of which 175 are 
at 0.05 < z < 0.3, where we expect NoSOCS to be very 
complete. 

The top panels of Figure [15] show the recovery rate 
of X-ray selected clusters by our survey as a function of 
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Fig. 15.— Top Left: NoSOCS recovery rate of NORAS clusters 
as f(zNORAs)- To P Right: NoSOCS recovery rate of NORAS 
clusters as f(Lx), for clusters with 0.05 < ZpfORAS < 0-3- Bot- 
tom Left: NORAS recovery rate of NoSOCS clusters as f(z p /iot). 
Bottom Right: NORAS recovery rate of NoSOCS clusters as 

{(Ngais)- 

NORAS spectroscopic redshift (left) and X-ray luminos- 
ity (right, for clusters with 0.05 < z NO ras < 0-3). At 
0.05 < z < 0.3 we recover 80-95% of the NORAS clus- 
ters, depending on the matching radius, with a distinct 
drop at z > 0.25, as expected from our completeness 
functions. At high Lx, we recover 100% of the NO- 
RAS clusters, but they are few in number. At moderate 
Lx (~ 10 45 erg s _1 ), the recovery rate is quite stable 
near 80%, mostly due to clusters missed at higher red- 
shifts. The bottom panels show the reverse compari- 
son, the recovery rate of optical clusters, as a function 
of NoSOCS photometric redshift (left) and optical rich- 
ness (right). As with other optical cluster surveys, we 
have nearly two orders of magnitude more candidates 
than NORAS, resulting in a very low recovery rate of 
NoSOCS clusters in the X-ray. The recovery rate in- 
creases with redshift, as the fraction of poor clusters de- 
creases in the optical. This effect is clearly illustrated in 
the bottom right panel, as the recovery rate approaches 
50% for clusters with N ga i s > 80. Nevertheless, it ap- 
pears that NORAS misses over half of the richest clusters 
in our sample. The extensive spectroscopy for NORAS 
also allows an additional check on our photometric red- 
shifts. The left panel of Figure [14] plots NoSOCS z phot 
against NORAS z spec , with matches within 0.5ft -1 Mpc 
shown in black, 0.5 < offset < 0.75ft _1 Mpc in green, and 
0.75 < offset < l.Oh^ 1 Mpc in red. For the 145 clusters 
matched within 0.5/i _1 Mpc, we find Qo-{zNoSOCS.,phot — 
ZNORAS,s P ec)/(l + z NO RAS,s P ec) = Az = 0.026, consis- 
tent with the errors estimated from the photometric re- 
lation combined in quadrature with the background es- 
timation errors. The recovery rates and typical o ffsets 
are also in good agreement with lLopes et al.l ([2006), who 
compared NoSOCS clusters to the more heterogeneous 
BAX database. 

6.2. X-ray Luminosities from RASS 

As seen in Figure 1141 X-ray measurements from NO- 
RAS are only available for a very small fraction of op- 
tically selected clusters. However, one can use the loca- 
tions of optically selected cluster candidates to measure 
X-ray fluxes and luminosities from RASS, allowing us to 
improve the optical-X-ray correlations. Even though the 
significance of X-ray emission in these areas may be too 
low to identify extended sources in RAS S, we can derive 
either fluxes or upper limits following Bohringc r et al.l 
(2000) to generate a much larger dataset. 

We first restrict the NoSOCS sample to clusters at 
0.069 < z p hot < 0.196, the same redshift range used for 
substructure measurement. The optical richnesses are 
most reliable at these distances, where the luminosity 
function is completely sampled over the magnitude range 
used to derive N ga i s . We further restrict the sample to 
clusters with N ga i s > 25, as we do for measuring R200, to 
avoid poor clusters/groups where the X-ray emission is 
unlikely to be detected or may be dominated by the X-ray 
halo of a single galaxy. The X-ray luminosities Lx are es- 
timated from count rates in ROSAT PSPC images taken 



as part of the RASS. Images and exposure maps in the 
0.4-2.4 keV band are retrieved from the ROSAT archive 
via FTP. We avoid the softest ROSAT channels (0.1- 
0.4 keV) since the background is highe r. To make our 
measu rements comparable to those of IBohringer et al.1 
(l2000h . we follow an almost identical procedure. The 
background is estimated in an annulus with an inner ra- 
dius of 20' and width of 21.3', divided into twelve sectors. 
The median count rate from these sectors is used as the 
background, after removing any sectors containing point 
sources . The details of this p rocedure are described in 
§3.1 of IBohringer et al.l (|2000h . The cluster X-ray flux 
is then computed using fixed apertures of 0.5 and 1.0 
hr 1 Mpc radius. We do not perform the growth curve 
analysis (GCA) because the vast majority of NoSOCS 
clusters have very low X-ray fluxes, making the GCA ex- 
tremely unstable. The computed X-ray fluxes are then 
corrected for flux missing from the fai nt outer regions 
using the technique described in §3.5 of [Bohringc r et aTl 
(2000). Finally, only clusters whose total counts are 3cr 
above the background are considered reliable. All oth- 
ers are reported as upper limits. The X-ray luminosities 
measured within fixed apertures of 0.5 and 1.0 h~ x Mpc 
(in units of 10 43 erg s _1 ) along with the associated errors 
and X-ray temperatures (in keV) can be found in Tabled 
Column 1 gives the cluster name, while Columns 2-4 give 
the X-ray luminosity, luminosity error, and temperature, 
all derived within a 0.5 h^ 1 Mpc aperture. Columns 5- 
7 provide the same quantities measured using a 1.0 h^ 1 
Mpc aperture. Luminosities marked with a ":" are upper 
limits. 

6.2.1. Validation with NORAS and REFLEX 

To test our methodology, we h ave recomputed Lx for 
the entire NORAS and REFLEX (|Bohringer et al.ll2001h 
samples using our software. We use the GCA-derived 
apertures reported in those two surveys, along with the 
redshifts, missing flux corrections and plasma models 
taken directly from the respective samples. Rather than 
transform between the different cosmologies used in NO- 
RAS and REFLEX, we perform all calculations with the 
cosmological parameters used in those surveys. To con- 
vert the measured total count rate to an unabsorbed 
X-ray flux in the full ROSAT soft energy band (0.1 
- 2.4 keV), we use the PIMMS tool available through 
NASA H EASARC. We assume a R aymond-Smith (RS) 
spectrum (Ra ymond fc Smith1ll977f ) to represent the hot 
plasma present in the intracluster medium, with a metal- 
licity of 0.2 of the solar value and the interstellar hydro- 
gen column density a long the line-of-sig h t tak en from 
iKalberla etldl (|2005l ) and iBajaja et all (j2005f l. The 
plasma temperature is estimated in two different ways. 
First, we use a fixed temperature of 5 keV, which is typi- 
cal for clusters (Markcvitch 1998), and term the resulting 
luminosity Lx5- Second, we use an iterative procedur e 
relying on the Lx — Tx relation from lMarkevitchl (|1998h . 
We start by calculating Lxs and finding the correspond- 
ing temperature, assuming an RS spectrum. This new 
temperature is used to recalculate the luminosity based 
on an RS spectrum, and the procedure is iterated un- 
til convergence is reached, when the change in temper- 
ature is ATx < IkeV, comparable to the scatter in the 
Lx — Tx relation. The procedure typically converges in 
two or three iterations. For both luminosity measures, 
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Fig. 16. — Comparison of X-ray luminosities computed using our 
methodology vs. NORAS (left) and REFLEX (right) . 



we apply a if-correction (|Bohringer et al.l [2000( 1 to de- 
rive the X-ray luminosity in the rest frame 0.1 - 2.4 keV 
band. 

The only significant methodological differences be- 
tween our technique and the previously published works 
are (i) we use the 0.4-2.4 keV images provided by the 
ROSAT archive, while they worked directly from the 
event files in the 0.5-2 keV range, (ii) we use a metal- 
licity of 0.2Zq instead of O.SZq, and (iii) they derive 
an independent count rate to flux conversion while we 
rely on PIMMS. Nevertheless, our results are in excel- 
lent agreement with both surveys. The comparisons to 
NORAS and REFLEX are shown in the top and bottom 
panels of Fig. [TBI respectively. We find very small offsets 
of ~ 5% in Lx between our measurements and the liter- 
ature values, likely due to differences in the count rate to 
flux conversion. The scatter is small, <r{5LogL x ) ~ 0.03 
over a very broad range of Lx, demonstrating that we 
are able to correctly recover X-ray luminosities with our 
technique. 

Four examples of the optical and X-ray properties are 
shown in Figure 1171 where we overlay contours of X-ray 
emission from the RASS on DPOSS .F-band images of 
four clusters in our catalog. The images are 0.5° on a 
side, centered on the NoSOCS optically-selected cluster 
centers. The clusters range in richness from N ga i = 30 — 
90, and redshifts of z = 0.07 - 0.183. There are 6 X-ray 
contours evenly spaced between the background level and 
2cr over the background. A circle of radius 0.5 x RAbeii is 
plotted, centered on the X-ray flux centroid. Although 
evident in the optical images, the X-ray fluxes are clearly 
not very high, and even moderately rich clusters near 
the median redshift of our catalog (such as the one at 
top right in the figure) do not stand out strongly. The 
X-ray contours are usually well matched to the optical 
center, except for the top left cluster. Visual inspection 
of the galaxy distribution in the latter field shows that 
the NoSOCS cluster center is in between two apparent 
overdensities which have been blended in our catalog, 
and only one of which is X-ray detected. This suggests 
that searching for clusters with highly discrepant optical 
and X-ray positions and/or fluxes can be used to find 
such projections. 

6.2.2. Optical vs. X-ray Properties 

The comparison of optical richness and our estimate 
of Lx from the iterative procedure described above is 
shown in Figure [TBI using 1649 clusters with 0.07 < z < 
0.19, Ng a i > 25 and successfully measured X-ray lumi- 
nosities with Lx > 5 x 10 42 ergs s _1 . Clusters where 
the X-ray luminosity is only an upper limit and those 
where the background estimation failed are not included. 




RA = 1*2 157623 DEC 37.816692 
Ngai = 63.25 zphol = 1831 



RA ~ 141.012039 DEC. 14.179153 
Ngai . 69.59 zphol . 0.1 372 



Fig. 17. — Contours of X-ray emission from the RASS overlaid 
on DPOSS .F-band images of four clusters in our catalog. The 
images are 0.5° on a side, centered on the NoSOCS optical cluster 
centers. There are 6 X-ray contours evenly spaced between the 
background level and 2a over the background. A circle of radius 
0.5 X RAbell is plotted at the X-ray flux centroid, while the optical 
center is marked with a cross. 

Individual clusters are plotted as dots, while the binned 
results (with each bin containing 200 clusters) along with 
their la scatter are shown as the large points with error 
bars. 

While the scatter is large , the binne d relat ionship 
agrees with that found by lLopes et al.l (|2006h using 
higher quality X-ray data, Lx cx N^ 16 . We also show, 
as the solid line, the relationship between Lx within 750 
kp c and A^no from the X-ray stacking analysis performed 
bv lRvkoff e t ail (|2008l) using the MaxBCG cluster cata- 
log (their equation 5) , but simply repl acing A^oo w i th our 
Nggi s- The powe r -law s lopes of the lLopes et al.l (2006) 
and iRvkoff et al.l (|2008[ ) relations are nearly identical, 
despite the completely different richness measures. As 
seen in Figure [THl a similar relationship holds for our 
clusters, despite ROSAT's limited spatial resolution and 
count rate as well as the limited photometric accuracy 
and depth of DPOSS, demonstrating that a reliable clus- 
ter sample can be defined from such data. It is also pos- 
sible transform our N ga i s to the MaxBCG A^on using 
Eqn. 7, and plot the relation from IRvkoff et all (|2008r ) 
using this pseudo-A^oo; this is shown as the dashed line 
in Figure [T8l We caution that this is not a reliable con- 
vers ion because the ric hness transformation is difficult 
and IRvkoff et al.l (|2008f ) use a completely different pre- 
scription for computing Lx- 

While the X-ray data from RASS is limited, especially 
for the poorer, lower mass systems, this catalog of in- 
dividual cluster X-ray measurements is the largest com- 
piled to date. It is only recently that astronomers have 
undertaken systematic comparisons of optical and X-ray 
cluster samples by returning to the source data and re- 
extracting physical properties consistent ly, rather than 
simply matching catalogs. For instance lDonahue et al.l 
(2002) compared independently detected X-ray and op- 
tical clusters from the same patches of sky. They found 
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TABLE 5 
X-ray Measurements 



Name 


Within 0.5/tff, 1 Mpc 




Within 1.0/ifn 1 Mpc 




Lix (,4U erg s ; 




T x (keV) 


L<x (Lu erg s ) 


FjTT(L,x ) 


T x (keV) 


nsc ti 1 1 7504-685Q1 n 


0.620 


0.140 


8 


0.160: 


0.260 


4 




0.200: 


0.960 


0.4 


1.160 


0.510 


1.1 


NSC 1091 71 1 -4-524442 


2.120 


0.600 


1.4 


28.760 


0.220 


5.3 


NSC 1102307+520201 


1.360 


0.830 


1.2 


5.020 


0.310 


2.2 


NSC 1101218+460643 


1.780 


0.520 


1.3 


1.020 


0.380 


1.0 


nsc n 038054-490496 


0.160 


0.030 


1 


0.160: 


0.620 


1 


NSC 11 7^1 54-37491 5 


0.360 


2.450 


0.6 


5.160 


1.320 


2.3 


NSC 1151120+363421 


1.590 


0.950 


1.2 


0.860 


0.910 


0.9 


NSC 1082043+301238 


0.120: 


0.000 


0.3 


1.050 


0.160 


1.0 


NSC 1152111+292632 


0.370 


1.920 


0.6 


2.900 


1.000 


1.7 


NSC 1081 Q494-9641 9Q 


0.020: 


0.090 


2 


3.190 


0.040 


1 8 


NSC 11 5531 94-973835 


0.570 


1.370 


7 


1.820 


1.370 


1 3 


NSC 109091 1 4-1 Q0446 


5.700 


0.250 


2.4 


8.510 


0.290 


2.9 


NSC 11 1 40474-1 81 Q39 


4.350 


0.580 


2 1 


2.060 


0.600 


l 4 


NSC J164837+193606 


0.270 


0.610 


0.5 


0.280: 


0.530 


0.5 


NSC J085246+161920 


0.800 


2.300 


0.9 


1.010: 


2.680 


1.0 


NSC J141229+140110 


4.480 


0.630 


2.1 


6.560 


0.520 


2.5 


NSC J011144+100349 


0.210 


0.120 


0.5 


1.700 


1.430 


1.3 


NSC J094338+085430 


0.190 


2.010 


0.1 


1.920 


1.540 


1.4 


NSC J135224+092048 


0.520 


0.030 


0.7 


4.010 


0.060 


2.0 


NSC J021010+080844 


0.690 


1.450 


0.8 


5.850 


1.370 


2.4 


NSC J104929+033846 


1.770 


0.800 


1.3 


2.530 


0.700 


1.6 


NSC J154555+030814 


0.930 


0.030 


1.0 


0.640: 


0.260 


0.8 


NSC J014426+021221 


0.570 


0.050 


0.7 


1.850 


0.040 


1.3 


NSC J104534-002506 


0.180: 


0.460 


0.1 


0.450: 


0.050 


0.7 


NSC J152156+013000 


0.240: 


0.000 


0.5 


1.870 


0.060 


1.1 



poor correlation between optical richness and X-ray lu- 
minosity, but could not pinpoint the physical reason for 
this, and pointed out the need to understand the effect of 
this sca tter on mass selectio n. The RASS-SDSS clusters 
survey (jPopesso et al.l2004l ) instead uses a small but very 
well measured sample of 114 X-ray detected clusters, and 
finds good correlation between X-ray luminosity or tem- 
perature and optical luminosity, if one has excellent data 
and chooses the measurement parameters (such as the 
aperture for richness measurement) carefully. However, 
it is worth noting that their sample rem ains one requir- 
ing X- ray detections, which was shown bv lDonahue et ail 
(2002) to potentially bias the results. 

The o nly other l arge o ptical - X-ray comparisons are 
those of iDai et al.1 (|2007[ ). who used stacking techniques 
to derive X-ray properties of over 4000 clusters selected 
opti cally from the Two-M icron All Sky Survey (2MASS), 
and iRvkoff et al.l (|2008( ), who stacked X-ray data for 
~ 17, 000 MaxBCG clusters from the SDSS. Even with 
the low redshift li mit (z < .1) im posed by the shallow 
depth of 2MASS, IDai et al.l (|2007f ) relied on stacking of 
X-ray data for clusters binned by their optical properties 
to measure correlations between mass (optical richness) , 
luminosity and temperature. They find similar corre- 
lations to those in the literature for individual clusters, 
but must model the Poisson fluctuations in the number of 
galaxies in a cluster of a given mass. At higher redshifts, 
where evolution in the cluster populations becomes more 
important, understanding and modeling these flu ctua- 
tions will be more challenging. IRvkoff et al. (2008) were 
able to examine some issues related to bias arising from 
scatter in the L^-richness relation with a small sample 
of clusters where individual X-ray measurements were 
possible. 

7. CONCLUSIONS 
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Fig. 18. — Comparison of optical richness N ga i and X-ray lu- 
minosity Lx within 0.5 h~ 1 Mpc for 1649 clusters with 0.07 < 
2 < 0.19, N ga{ > 25 and L x > 5 X 10 42 ergs s _1 . Individual 
clusters arc plotted as dots, while the binned results (with each 
bin containing 200 clusters) along with their rms errors shown as 
the large squares with error bars. The solid line shows the relation 
Lx,7S0k P c = e 3A (§ found by IRvkoff et all l|200ST ). directly 
replacing their N200 with our N ga [ s . The dashed line shows the 
same relation, but now transforming our N ga [ a to MaxBCG N200 
using Eqn. 7. Dotted lines show the N ga i > 25 limit for attempt- 
ing to measure Lx and the Lx > 5 X 10 42 ergs s" 1 limit imposed 
on the sample when computing binned medians. 

We have presented NoSOCS, a new cluster catalog 
based on the \b\ > 30° plate scans from the Digi- 
tized Second Palomar Observatory Survey. Spanning 
over 7T steradians, this is the largest area o ptical clus- 
ter catalog created sinc e those of lAbelll Jl958) and 
lAbell. Corwin fc Olowinl (1989). In terms of area cover- 
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age, it will only b e superseded by new sky surveys suc h 
as Pan-STARRS (|Kaiserll2004ft and LSST (jTvsonll2006r ). 
both of which have cosmology through clusters as im- 
portant science drivers. We show consistency among the 
three regions covered by NoSOCS, and with th e SDSS 
MaxBCG cluster catalog of lKoester et all ( |2007bh . How- 
ever, interesting discrepancies between these two large 
surveys remain. These include large numbers of poor 
clusters missed by one survey but found in another, 
suggesting lower completeness, higher contamination, or 
some combination of the two for such systems, in cither 
or both surveys. Even for supposedly rich clusters there 
are sufficient discrepancies to call into question our abil- 
ity to use such surveys for high-precision cosmological 
constraints. Understanding the sources of these disagree- 
ments requires further investigation of both systematic 
errors and individual cluster candidates. 

We have also derived X-ray luminosities for a large 
subset of our cluster sample from ROSAT all-sky X-ray 
survey data. We demonstrate that the optical richness 
and X-ray luminosities are well correlated, albeit with 
moderate scatter. We find that, despite the poor photo- 
metric data and low X-ray luminosities of most NoSOCS 
clusters, the correlation between L x and N ga i s is in good 
agreement with literature results using better data and 
stacking analyses. Refinements to both the optical rich- 
nesses and especially deeper X-ray survey data will be 
necessary to improve this relation and truly understand 
the utility of optical richnesses for mass estimation. Nev- 
ertheless, our results show promise for using large surveys 



for such measurements in cases where the data quality 
is less than superb, as may be expected for the high- 
est redshift clusters even in upcoming deep surveys such 
as Pan-STARRS and LSST. Furthermore, our ability 
to measure X-ray luminosities for hundreds of clusters 
not originally detected in the RASS argues for improved 
multi- wavelength detection methods that leverage multi- 
ple surveys (optical, infrared, X-ray, S-Z) to find distant 
and/or poor clusters which would otherwise fall below 
the significance cutoff in a single passband. 
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